Fractional impurity moments in two-dimensional non-collinear magnets 
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We study dilute magnetic impurities and vacancies in two-dimensional frustrated magnets with 
non-collinear order. Taking the triangular-lattice Heisenberg model as an example, we use quasi- 
classical methods to determine the impurity contributions to the magnetization and susceptibility. 
Most importantly, each impurity moment is not quantized, but receives non-universal screening cor- 
rections due to local relief of frustration. At finite temperatures, where bulk long-range order is 
absent, this implies an impurity-induced magnetic response of Curie form, with a prefactor corre- 
sponding to a fractional moment per impurity. We also discuss the behavior in an applied magnetic 
field, where we find a singular linear-response limit for overcompensated impurities, and propose 
experiments to test our theory. 



Impurities have been established as a powerful means 
to both probe and tune bulk properties of correlatcd- 
electron materials. In quantum magnets, non-trivial phe- 
nomena include vacancy-induced magnetism in quantum 
paramagnets [1] and quantum percolation [2]. Single- 
impurity behavior has been predicted to be exotic in 
quantum critical magnets, where a universal fractional 
Curie moment appears at low temperatures [3-5]. Iso- 
lated impurities in magnets with long-range order have 
been studied as well, with most works focussing on the 
square-lattice Heisenberg magnet [3, 4, 6, 7]. 

This paper is devoted to impurities in geometrically 
frustrated spin-5 magnets which order non-collinearly - 
a topic which has received little attention [8] . As we show 
below, vacancies (i.e. non-magnetic impurities) in non- 
collinear magnets display a behavior which is richer and 
qualitatively different compared to their collinear coun- 
terparts. In particular, the magnetic moment m associ- 
ated with a single vacancy is not quantized, in contrast to 
the collinear case [3] where it is locked to m = S. This ef- 
fect is already present at the classical level: nearby spins 
re-adjust their directions in response to the vacancy, re- 
flecting that frustration is locally reduced. This partially 
screens the vacancy moment, with the screening cloud 
decaying algebraically due to Goldstone modes. 

At zero temperature, the direction of the vacancy mo- 
ment m is fixed by the bulk magnetic order. In contrast, 
at T > in two dimensions (2d) there is no long-range or- 
der due to the Mermin- Wagner theorem, and the vacancy 
moment is free to rotate. This rotation is classical, as it 
is coupled to a rotation of the bulk spins surrounding the 
vacancy [3, 4]. As a result, the linear-response suscepti- 
bility has a singular piece, Ximp(T) = m 2 /(3fcT), corre- 
sponding to the Curie response of a fractional moment 
for each vacancy [9]. For the triangular-lattice Heisen- 
berg antiferromagnet (AF) with nearest-neighbor inter- 
actions, we find in a 1/5* expansion 



described in detail below. In stark contrast to fractional 
effective moments found at bulk or boundary quantum 
critical points [3, 5, 10], the present mechanism is realized 
deep inside the renormalized classical regime [11] of a 2d 
magnet. 

A finite magnetic field h has two effects which tend 
to compete: it orients the impurity moment parallel to 
the field and it induces a macroscopic bulk moment. 
This bulk-boundary competition is governed by a field- 
induced length scale lu oc 1/h and limits the linear- 
response regime [12]. For a single overcompensated im- 
purity in the triangular lattice, we find this competition 
to be particularly drastic: Linear response breaks down 
at any finite field. 

In the body of the paper, we sketch the derivation 
of these results and propose tests and extensions of the 
non-trivial screening advocated here. Our considerations 
qualitatively apply to a large class of frustrated AFs 
with non-collinear ground states, which are unique up 
to global spin rotations. For definiteness, we will present 
results for the spin-S* triangular-lattice Heisenberg model 



(2) 



m = -0.0405 + 0.196 + 0(1/ S) , 



(1) 



where a negative sign corresponds to overcompensation, 



The biquadratic exchange, with its strength parameter- 
ized by k = K/(JS 2 ), generates a family of models and, 
in particular, lifts the accidental classical degeneracy of 
the simple nearest-neighbor Heisenberg model in an ap- 
plied field [13]. In zero field, the ground state is given 
by the familiar coplanar 120° ordering at wavevector 
Q = (4vr/3, 0) for -2/9 < k < 2/9. 

Vacancy in the ground state of a classical non-collinear 
magnet. Consider a bulk AF with geometric frustration, 
where not all energetic constraints (e.g. all neighboring 
spins pairwise antiparallel) can be satisfied. Removing 
a single spin locally reduces frustration due to the elim- 
ination of constraints. For a non-collinear magnet, this 
seeds a re-adjustment of spin directions. 

For the triangular lattice, this re-adjustment is illus- 
trated in Fig. la, which shows the result of a numerical 
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FIG. 1: a) Classical ground-state spin configuration of the 
triangular AF with vacancy, showing the spin re-adjustment 
near the vacancy; the dashed arrows indicate the original 120° 
order, b) Rotation angle |<50|(r) along a high-symmetry line 
for k — and L — 51, together with the asymptotic power law 
SO oc l/r' ! . The inset shows the same data (solid), together 
with \S@(r)\ at finite field h/J = 0.5 (dashed), 1.0 (short 
dash), 1.5 (dash-dot), all for k = —0.05, for one sublattice 
in a log-linear plot. The exponential (instead of power-law) 
decay for h > is obvious. 
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FIG. 2: a) Effective vacancy moment m/S for the classi- 
cal triangular Heisenberg AF as function of the biquadratic 
exchange k in Eq. (2). b) MC results for Xim P as func- 
tion of T/S 2 , calculated with JS 2 = 1 and k = 0. The 
dashed line shows the predicted Curie law m 2 /(3kT) (3) with 
|m/S| = 0.04. The inset shows the low-temperature data in a 
log-log plot. The data in a (b) have been obtained for systems 
of size L = 51 (L — 9, 12); finite-size effects are negligible. 



optimization of the spin directions in a finite system of 
size L 2 with L = 51, where a single spin has been re- 
moved. The spins remain coplanar and rotate by angles 
50 relative to the original 120° configuration, such that 
the spins near the vacancy tend to be more antiparallcl. 
The numerical results for SQ show a sixfold (/-wave) an- 
gular symmetry and are consistent with a spatial decay 
of 6Q(r) oc 1/r 3 , Fig. lb [8]. 

This result is rationalized as follows: The full rotation 
pattern can be understood as the response of the sys- 
tem to a field h which couples to the six neighbors of 
the vacancy such that these spins are rotated towards 
an antiparallcl configuration. Hence, the field h acting 
on these six sites is locally transverse and alternating, 
in a rotated frame compactly written as h^lj=iPjSj 
with (3j = (— iy. The long-distance rotation is deter- 
mined by a transverse susceptibility, which is dominated 
by the modes near the ordering wavevector with lin- 
ear dispersion uj q . A straightforward calculation gives 
5Q(r) cx / d d qe^ p l3 q /uj 2 cx l/r d+1 . 

The state with a single vacancy has a finite magne- 
tization to. While this would simply be to = 5 with- 
out re-adjusted angles (i.e. in the collinear case), the 
re-adjustment tends to screen this moment. For the tri- 
angular lattice, the numerical result, obtained from inte- 
gration over the screening cloud, is m/S — —0.0396(3), 
i.e., the missing spin is cwercompensated, such that the 
total moment points in the direction of the removed spin. 
The value of to is non- universal, i.e., depends on details 
of the Hamiltonian: Fig. 2a shows m/S as function of 
the biquadratic exchange coupling K in Eq. (2). 

Vacancy: 1/S corrections. Quantum corrections to the 
classical T = results can be obtained using spin-wave 
theory. Holstcin-Primakoff bosons a are introduced to 
capture deviations from the classical state in the presence 



of a vacancy, Fig. la. Upon expressing the Heisenberg 
model in terms of the a bosons, terms linear in a vanish 
as required. Linear spin-wave theory amounts to a diago- 
nalization of the quadratic-in-a piece of the Hamiltonian, 
which has to be done numerically for finite lattices [14] 
due to the inhomogencous reference state. Doing so, we 
find the full spectrum of eigenvalues and eigenvectors, 
which can be used to calculate 1/5 corrections to ther- 
modynamic observables as well as response functions. 

The local magnetization correction 5m(fi) — (ajcii) 
decays to the known bulk value of Smt = 0.26 [15] at 
long distances, corresponding to a staggered magneti- 
zation of m\> = 5 — 0.26. The impurity contribution, 
8m{fi) — Smii, indicates enhanced quantum corrections 
near the impurity which fall off as 1 /r 3 , consistent with 
the Goldstone-mode expectation. The 1/5 correction 
to the uniform moment associated with the vacancy is 
obtained from integration, 5m = Sm{fi) cos Q(f f i), 
which evaluates to 6m = 0.196, Eq. (1). Further correc- 
tions at higher orders in 1/5 will not qualitatively modify 
the result of a non-universal fractional value of to, but 
apparently both overcompensation and undercompensa- 
tion may occur, depending on 5 and microscopic details. 
We note that local impurity-induced magnetization cor- 
rections obeying Sm(fi) — Smi, oc 1/r 3 also occur in the 
collinear square-lattice case, but here spin conservation 
demands that the integral 8m vanishes; hence to remains 
locked to 5 [4, 6]. 

Finite temperatures: Fractional Curie response. For 
T > in two space dimensions, long-range bulk magnetic 
order is destroyed by thermal fluctuations, with the cor- 
relation length £ being exponentially large at low tem- 
peratures, T <C J. Consequently, the direction of the 
impurity moment is no longer fixed, but is free to rotate 
with the local orientation of the bulk magnetic domain 
surrounding the impurity. It has been shown both an- 
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alytically [3, 4] and numerically [7] that this rotation is 
classical and leads to a linear response of Curie form: 



a) y 



Ximp(T) 



3fcT 



+ O(T ) , 



(3) 



where the subleading term receives a multiplicative log- 
arithmic correction in 2d [4, 7, 16]. 

In the non-collinear case, the partial screening of the 
vacancy moment, established above for T = 0, will re- 
main intact at small T > because of the large corre- 
lation length. This implies a Curie response (3) corre- 
sponding to a fractional moment per vacancy. This cen- 
tral result is fully borne out by numerics: we have per- 
formed classical Monte Carlo simulations of triangular- 
lattice Heisenberg magnets, using the standard Metropo- 
lis algorithm. In Fig. 2b we show the result for the im- 
purity susceptibility Ximp, obtained from subtracting the 
linear-response % of a system with vacancy from that of 
a system without vacancy. While the high-temperature 
part is difficult to analyze given the error bars, the low- 
temperature data clearly show a Curie divergence, with a 
prcfactor consistent with |m/5| = 0.04 within error bars. 
For S < oo, we expect a subleading logT contribution 
to Ximp(T) arising from Goldstone modes similar to the 
collinear detailed analysis will appear elsewhere. 

Vacancy vs. extra spin. So far, we have considered 
the special case of a vacancy, experimentally obtained by 
replacing a magnetic by a non-magnetic ion. A different 
type of impurity is an extra spin of size S' , coupled to a 
single site of the bulk magnet with a Heisenberg coupling 
J'. For antifcrromagnctic J' ^> J and S — S' , the impu- 
rity spin and its bulk partner lock into a singlet, and we 
recover the vacancy case. On the other hand, for J'c J 
the re-adjustment of the spin directions due to the impu- 
rity will be parametrically small in J'/ J, and we expect 
for the impurity moment m — ¥ S as J' — >■ 0. Hence, vary- 
ing J'/ J leads to a continuous change of m. Interestingly, 
this crossover cannot be captured in a 1/iS expansion for 
the extra-spin problem: In the classical limit, there is no 
singlet formation for large J', and this is not recovered 
at any order in 1/5, as can be seen by explicit calcu- 
lation. However, the direct vacancy calculation can be 
performed in l/S, as shown above, the crucial difference 
being that the missing spin (or the J' = oo singlet) is 
built in from the outset. 

Finite magnetic field. For the collinear square-lattice 
AF, it has been shown that a vacancy in a finite applied 
field generates spin textures in its vicinity, which result 
from the competition between aligning the vacancy mo- 
ment and inducing a bulk moment [12]. Here we investi- 
gate the non-collinear case on the triangular lattice. As 
the nearest-neighbor Heisenberg model has an acciden- 
tal degeneracy of classical ground states at finite fields, 
which is lifted in favor of coplanar states (Fig. 3a) both 
by quantum and thermal fluctuations [17, 18], we choose 
to investigate the classical model with biquadratic ex- 
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FIG. 3: a) Schematic evolution of the three-sublattice copla- 
nar bulk spin configurations as function of applied field for the 
triangular AF, with the 1/3 magnetization plateau shaded. 
These states are selected out of the classical k = ground- 
state manifold of H (2) by both thermal and quantum fluctu- 
ations as well as negative small k. A single vacancy chooses 
one of the sublattices: b) undercompensated and c) overcom- 
pensated case. In b), an impurity quantum phase transition 
occurs inside the plateau phase. (Note that undercompen- 
sation does not occur in our family of classical models with 
k < 0, but is expected for S < oo from Eq. (1).) d) Spin 
configuration with (overcompensated) vacancy, calculated for 
h/J = 1.0, k = -0.05, and L = 51. 



change, Eq. (2) with —2/9 < k < 0, which leads to the 
same coplanar finite-field phases as the ones selected by 
fluctuation effects. (Non-coplanar states are favored for 
< k < 2/9.) 

Qualitatively, the vacancy physics strongly differs 
between the undercompensated and overcompensated 
cases. For undercompensation, Fig. 3b, a small field will 
orient the system such that the vacancy sublattice points 
antiparallel to the field. This is compatible with the field- 
induced bulk state, hence a strong competition between 
bulk and boundary effects is absent, and the zero-field 
limit will be smooth. 

This is different in the overcompensated case, where 
orienting the vacancy moment in field direction is incom- 
patible with the bulk state. Our numerics shows that the 
system chooses a compromise such that the vacancy sits 
in one of the sublattices directed approximately parallel 
to the field, with a significant distortion near the vacancy, 
Fig. 3c, d. This distortion falls off exponentially (there is 
no coupling to the remaining Goldstone mode), with a 
length scale lh <x 1/ft [12], Fig. 4b. Most importantly, the 
zero-field limit is singular in this case, i.e., the distortion 
pattern for ft — ^0 does not recover its zero-field structure. 
This is seen in both the inset of Fig. lb and Fig. 4a. The 
latter shows mj mp (ft), defined as the difference of the to- 
tal magnetizations with and without vacancy. By con- 
struction, mi mp (ft = 0) = \m\ and mi mp (ft — > oo) = —S. 
Fig. 4a demonstrates that mi mp (ft— >0) again represents 
a fractional impurity moment which is different from |m| 
in the overcompensated case. 

The evolution of mi mp (ft) through the bulk magneti- 
zation plateau is also very different in the undercompen- 
sated and overcompensated cases, Fig. 3b, c. While it is 
smooth for undercompensation, a jump occurs at ft = 3 J 
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FIG. 4: a) Impurity contribution to the magnetization, 
w,i mp (/i)/S, as function of applied field h, for the classical 
triangular AF (2) with k = —0.05. The shaded region cor- 
responds to the bulk 1/3 magnetization plateau. The inset 
shows a zoom onto the small-field region; squares (crosses) 
are data for L = 51 (L = 21). The circle at h = rep- 
resents the linear-response value \m/S\, demonstrating the 
breakdown of linear response for this overcompensated case, 
b) Field-induced length scale lh(h) obtained from an expo- 
nential fit to 8Q(r) for L = 51, together with the anticipated 
lh on 1/h behavior (dashed). Finite-size effects are important 
for small h > where l h <C L is violated. 



in the overcompensated case, Fig. 4a. This signals a first- 
order impurity transition where the vacancy site switches 
the sublattice implying that the entire bulk configuration 
rearranges. 

Finite impurity concentration. We finally discuss the 
measurable consequences of our findings in the realistic 
case of a finite impurity concentration 7ij mp . Assuming 
that the impurities are distributed equally over all sub- 
lattices, their moments tend to average out at T = and 
h = 0. This behavior persists at finite T, provided that 
n imp i s the mean impurity dis- 



£ > limp where 

tance. In the opposite limit, £ <C l lmp , the impurity mo- 
ments fluctuate independently, and their response simply 
adds up. Hence, observing fractional Curie response of 
independent impurity moments is possible at elevated T 
and small ni mp [9] . Note that elevated fields which induce 
hi *C ^imp also lead to an effective decoupling of multi- 
ple impurity moments, which, however, are polarized in 
this limit. The spin re- arrangement predicted to occur 
inside the plateau phase for overcompensated impurities 
is detectable by local probes like NMR. 

Conclusions. For impurities in non-collincar magnets, 
our main result is a partial screening of the impurity mag- 
netic moment, leading to a fractional Curie response at 
low temperatures in the 2d case. We have evaluated the 
vacancy moment for the spin-S 1 triangular-lattice AF in 
a 1/S expansion, but we expect our qualitative results to 
be valid for any frustrated AF with non-collinear ground 
state (which is unique up to global spin rotations). 

Our predictions could in principle be verified by large- 
scale numerical studies in analogy to Refs. [5, 7], how- 
ever, quantum Monte-Carlo approaches are plagued by 
the sign problem which is serious for most frustrated AF. 



On the experimental side, one can expect the physics 
described here to be generically realized, as Curie tails 
in x{T) due to impurities arc routinely observed in mag- 
nets. A quantitative analysis of these tails in samples 
with known concentration of impurities would allow to 
extract the fractional moment size m (in a regime where 
interactions between the impurity moments are small); 
our prediction is m C 5 in contrast to the behavior in 
collinear magnets. 

An interesting open question is how the fractional mo- 
ment advocated here evolves upon approaching a quan- 
tum critical point of the bulk magnet, where at criti- 
cality a universal fractional response is expected. Our 
results also call for investigations of vacancies in frus- 
trated collinear magnets where vacancies may induce 
non-collinear spin textures in order to reduce frustration. 
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